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Circular  Array  of  Dipoles 
Above  a  Perfectly  Conducting  Cylinder 


l.  INTRODUCTION 

Recent  developments  in  low-inertia  phase  shifters  and  electronic  switches^ 
have  aroused  Interest  in  large,  circular  arrays  on  conducting  cylindrical  surfaces 
for  electronic  agile  beam  positioning  or  radar  resource  optimization,  invariant 
beam  characteristics  over  the  entire  360°  azimuthal  scan  sector,  and  frequency- 
independent  beam  pointing  direction.  A  low-sidelobe,  high-gain  aperture  design 
for  a  circular  array  of  dipoles  surrounding  a  perfectly  conducting  cylinder  requires 
a  knowledge  and  control  of  mutual  coupling  parameters  in  the  array  environment.^ 

3 

Carter  in  1943  outlined  a  rigorous  solution  of  Maxwells'  equations  for  a  dipole 
near  a  long  cylinder  (see  Figure  1),  and  tabulated  formulas  for  the  radiation  pat¬ 
terns  for  three  different  circular  array  configurations.  Carter  used  the  reciprocity 
theorem  and  the  k  >own  solution  for  scattering  from  a  cylinder  to  obtain  far-zone 

(Received  for  publication  30  September  1980) 

1.  Provencher,  J.H.  (I9f0)  A  survey  of  circular  symmetric  arrays.  Phased 

Array  Antennas  (Oliver  and  Knittel,  Artech  House,  Inc. ,  Dedham,  Mass. ) 
page  293. 

2.  Hessel,  A.  (1970)  Mutual  coupling  effects  in  circular  arrays  on  cylindrical 

surfacea-«perture  design  implications  and  analysis.  Phased  Array  Antennas. 
(Oliver  and  Knittel,  Artech  House,  Inc,,  Dedham,  Mass. )  page  3y3. 

3.  Carter,  P. S.  (1943)  Antenna  arrays  around  cylinders.  Proceedings  of  IRE, 

Vol.  31,  December, 
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Figure  1.  Dipole  Near  a  Long  Cylinder 
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patterns.  Lucke  '  ‘  used  Green's  function  method  to  yield  expressions  for  the 
fields  in  terms  of  an  integration  in  the  complex  plane.  The  problem  was  also  solved 
by  Harrington'*'  '  by  converting  a  three-dimensional  radiation  pv-oblem  having  cy¬ 
lindrical  boundaries  into  a  two-dimensional  problem  by  applying  a  Fourier  trans- 
form  with  respect  to  the  cylinder  axis,  Z.  Neif  also  used  the  Fourier  transform 
method,  and  by  means  of  a  limiting  process,  obtained  an  approximate  form  for  a 
dipole  near  a  finite  cylinder. 

In  this  report,  the  problem  of  finding  the  electromagnetic  field  and  admittance 
of  a  circular  array  of  dipoles  near  and  parallel  to  an  infinite  conducting  cylinder 


4.  Lucke,  W.  (1949)  Electric  Dipoles  in  the  Presence  of  Elliptical  and  Circular 
Cylinders,  Heport  No.  1,  Project  188,  Stanford  Research  Institute, 

Ktanford,  Calif. 

r>,  Lucke,  \V.  (1951)  Electric  dipoles  in  the  presence  of  elliptical  and  circular 
cylinders,  Joui  nal  of  Applied  Physics,  Vol.  22,  No.  1,  January. 

t).  Harrington,  H.  F.  .and  l.ePage,  W.  R.  (1949)  A  Study  of  Directional  Antennas  for 
DIF  Purposes,  Report  1,  Departmeni  of  Electrical  Engineering,  Syracuse 
I’niversity,  Syracuse,  NY,  September. 

7.  Harrington,  R.  ' (I9(il)  Time-Harmonic  Electromagnetic  Fields.  McGraw-Hill, 

New  Y’ork,  N")',  Chapter  3. 

8.  Neff,  H.  P. ,  Hickman,  C.  E.,  and  Tillman.  J.D.  (1964)  Circular  Arrays  Around 

Cylinders.  Report  No,  7,  Department  of  Electrical  Engineering,  University  of 
Tennessee,  Contract  AFl9(628)-288,  June. 
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has  baen  mathematically  formulated  Into  Integral  equatlona  via  the  Fourier  trane* 
form  method,  which  relates  the  Illumination  function's  voltages  to  the  resultant 
current  distribution  due  to  the  mutual  coupling.  We  use  the  computer  to  numerically 
solve  the  m  nonhomogeneous  Integral  equations,  using  numerical  approximations 
for  the  Integrals  and  complex  matrix  techniques  to  solve  the  simultaneous  equations. 
Such  an  analysis  should  account  for  both  the  space  wave  with  its  grating  lobe  and 
the  creeping  wave  resonances  on  the  curved  surface,  and  have  an  impact  on  the 
development  of  optimum  methods  tor  determining  techniques  for  radiation  pattern 
generation  with  very  low  sldelobes,  optimum  spacing  and  a  minimum  number  of 
array  elements  isee  Figure  2).  We  snail  Investigate  the  effects  of  mutual  coupling 
for  a  circular  array  located  a  quarter  of  a  wavelength  above  a  perfectly  conducting 
cylinder  having  a  radius  of  two  wavelengths.  This  simple  case  was  chosen  to 
easily  investigate  the  m  almuUaneotis  Integral  eo,uatlons  and  to  lay  the  foundation 
for  fliture  Investigutions  into  the  more  practical  case  of  a  perfectly  conducting 
cylinder  with  radius  of  approximately  10  wavelengths. 


OeaECTiVE  ,  TO  OBTAIN  SBO-OEGREE  SCAN  CAPABILITY  FOR 
LOW  SIOELOBE  ARRAYS  OF  ELEMENTS 


I.  MUTUAL  COUPLING  EFFECTS 

Z.0PTIMUM  ILLUMINATION 
FUNCTION 


I  LOW  SIOELOBE  AZIMUTHAL 
PATTERN 

2.CREEPING  WAVE  PHENOMENA 


I.  ELEVATION  PATTERN  CONTROL 

^STACKIN6  OF  CIRCULAR 
ARRAYS 


APPLICATIONS.  UNATTENDED  RADAR,  GROUND  BASED  TACTICAL  RADAR 


Figure  2.  Circular  Arrays  Around  Cylinder,  Optimization  Approach 


2.  MATHEMATICAL  ANALYSIS 


The  clrculcr  antenna  array  to  be  considered  (see  Figure  3)  consists  of  identi¬ 
cal,  parallel,  cylindrical  dipoles  equally  spaced  around  the  circumference  of  a 
circle,  which  is  concentric  to  the  perfectly  conducting  cylinder  which  the  array  of 
dipoles  surrounda.  Only  center-fed  dipoles  with  a  half-length  X/4  and  a  radius  b 
will  be  used  as  elements.  The  axial  center  of  each  dipole  is  perpendicular  to  the 
plane  of  the  circle.  Each  of  the  elements  of  the  array  is  thus  in  the  same  geomet¬ 
rical  environment.  The  radius  of  the  perfectly  conducting  cylinder  Is  a.  and  the 
radius  of  the  circular  array  of  dipoles  is  r^  (see  Figure  4).  The  expression  ior 
the  vector  potential  A  (r,  4,  s)  has  been  developed  for  a  circular  array  of  m 

Q 

dipoles  located  about  a  perfectly  conducting  cylinder. 

For  a  current  density  j  located  above  a  perfectly  conducting  cylinder  (see 
Figure  5), 

V*  A  +  A  =  .  (1) 

For  a  unit  dipole  at  (r  .  4  ,  z  ), 
o  o  o 

-  .  6(r  -  r  )6(4  -  e  )6(z  -  z  ) 

j  =  z  - 2 - _2 - 2.  .  (2) 

Therefore,  we  need  only  be  concerned  with  the  component: 

00 

f(z)  F(o.  r.  4)  do  (3) 

-00 

where 

00 

F(o.  r.4)  =  J]  .  (4) 

I'A  =  -00 


9.  Fante,  Dr.  Ronald  L.  (private  communication). 
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cyuMomCAL  boundary 
CONTAINING  CURRENT 
CLEMENTS 


Figure  4.  Perspective  View  of  Dipole,  Conductive 
Cylinder  Geometry 
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Figure  5.  Geometry  for  the  Development  of  Integral  Equations  for  a  Circular 
Array  Around  a  Perfectly  Conducting  Cylinder 


From  Eqs.  (3)  and  (4)  we  get 


In  cylindrical  coordinates,  the  laplacian  is 


Substitute  Eqs,  (f>),  (5),  and  (2)  into  Eq.  (1)  to  yield 


ir  E 


,i<»z  y  I  3 


r  SF  I  9r 


^  .2  2  m-  \ 

+  K  -  Or - K—  1  a 


6(r  -  r„ )6(*  -  4„)6(z  -  e  ) 


s 


12 


®n  ■  ‘^n 

n  n 


-ftjjfer  • 


Substitute  Eq.  (IS)  into  Eq.  (14)  to  yield 


Define  ^  =  J_(x)  H*(y)  -  J'(*)  H„(x)  (Wronskian),  then 
n  n  n  n 


-m„H„(x) 
o  n 

2»  X  A 


Also  from  Eq.  (IS), 


(15) 


(10) 


-M  J„(x) 
B  -  n  •  Q 
"n  “n  2xX  A 


(17) 


Now  we  need  to  satisfy  the  boundary  conditions  at  the  surface  r  -  a  of  the  cylinder, 
which  is  that  the  electric  field  tangential  on  its  cylindrical  surface  must  be  zero. 


E  = 


vx  vx  A 


,  9A\  ,  3A, 

37^ 


VX  VX A  =  r 


a^A  .  ,  3^ A, 

_ i  a.  I  1  *  j.  S 

drdz  *  r  9zd4 


1  3 

.  1 

3^A 

z 

■  r  5r 

V  / 

T 

3^^ 

r 

therefore 


dzd4 


=  0 


r=a 


(18) 


13 


14 


“n  -  ~H;;!y) 


Upon  8tibstitutin£  Eq.  (16)  Into  Eq4  (24)  and  rearranging  terms  we  get 


>11  ■  iixa  j 


Using  Eq.  (25)  in  Eq.  (17)  we  get 


B,  •  SoT^T?! 


Since  »  X  A  =  2i  we  get 


Cn  -  - 5 


D 


n 


M„H„(x) 
*^0  n 

4i 


Jn<y> 

Hn(y) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


where 


r^  =  radius  at  which  dipole  is  located 
a  =  radius  of  cylinder. 
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fo'"  r  2  we  i» 


E  f  -o  (a  /  t 

in=l  -L/2  -«  n*~« 


X  r)  =  1  co«  3^  +  ^.In  p\z\  ]  (33) 

where  the  dipole  la  at  (r,  d.  e).  Thus,  for  an  m-dlpole  circular  array,  we  have 
in  slmultaneoua  equations  to  solve.  Let  a  -  */23,  then  Eq.  (33)  becomes 


m  L/2  oe 

AE  /  ■’•oW/ 

m  =  l  -L/2  -« 


^lo(e/23.z^)  ^  ^ 
n=-*e 


B„(a)H_(‘^ r)  = 


n  n 


(34) 


»2 


Substitute  Eq.  (27)  into  Eq.  (34)  and  let  1  (a  )  =  A  cos  -r^  then, 

tn  o  ni  L 


m 


L/2 

A  E*m  / 

m^l  -L/2 


dz  cos 

o 


*o  f 

L  J 


do  e 


n-  -go 


r  J  (x)H  (y)  -  H  (X)J  (y)  /-s - 5  T  -4»V-. 

N/  I  _2 - EL- _ □ _ BJ_  ..  i»/..2  2  I  K 

""L  H(y)  • 


(35) 


Since  all  dipoles  are  equidistant  from  the  perfectly  conducting  cylinder,  r  -  r 


then, 


-  0“  r)  II  (V^  -  J  r  )  =  H 


o  n 


(X) 
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which  from  symmetry  becomes 


m 


E 

m  =  l 


^ln(*^+b/r^-4^)  j? 

0 


cos 


4 


X 


J_(x>  H  (x) 
n  n 


(37) 


11  Mf  asi  MiNimI  AaMittssc* 

Now,  let 


oe 

Tkm<“>  =  51  ^ 

n=0 


m(4^+b/r^-^^)  J  da  ^ 


Lo 

X_  co^  ^ 

o  4 


X 


J„(x)H.(x) 
n  Tt 


H„(,) 


then 


nt 

E 

m  =  i 


-j£t 


Now  we  shall  consider  the  M  iitesral  equations  by  taking;  a  look  at  the  equations 
yielded  for  each  base  voltage  Vj^  (K  -  I,  ...  m)  simultaneouBly.  We  get 


AjT^l  +  A2TJ2  *  "•  ■^m^lm  '  iTc 
A,T„  ...  A„Tj„  .  ^  Vj 


A,T,„i  +  A,T^,  +  . . .  A,„T  „  -  — ^ 

1  ml  2  m2  111  mm  u  ^  m 


(38) 
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tn  a  matrix  notation,  Eq.  (38)  ia  equal  to 


Til  "^12  ••• 

T  "^22  *  ‘  * 


'’'ml  ^m2  “■ 


• 

*  - 

• 

■’'im 

^2m 

a 

a 

s 

a 

a 

T 

mm 

a 

A 

m 

^m 

where 


(39) 


Therefore, 

fXV  . 


Given  valuea  for  Vj,  Vj,  ....  and  having  numerically  approximated  the  inte¬ 
grals  (m  =  I . m;  K  =  I . m),  we  compute  the  values  for  the  A^^'s 

by  factoring  the  matrix  T  into  the  L  -  U  decomposition  of  a  row -wise  permutation 

of  T  and  solving  the  systems  of  Eq.  (37).  Once  the  values  for  the  A  's  have  been 

m 

found,  we  can  then  determine  the  mutual  and  self  admittances  of  the  circular  array. 
Equation  (33)  gives  the  relationship  between  the  currents  and  voltages  of  the  array 
elements.  For  an  m-element  array. 


(40) 


Ym^l  ^  ^^2^2 


^  ^mkVk  ^ 


+  V  =  I 
mm  m  m 


where  Yu,  ¥22,  •  •  •  Yj^j^,  . . .  Y^^^  are  the  self  admittance  of  the  respective  ele¬ 
ment,  is  the  mutual  admittance  between  a  and  b,  is  the  applied  voltage  on 
the  element,  and  is  the  current  in  the  element. 

For  an  array  with  a  single  active  element,  which  for  simplicity's  sake  we  let 
be  element  no.  1,  every  Vj^  where  K  =  1  is  equal  to  zero.  Thus,  Eq.  (16)  becomes 
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h  •  ^11^1 

*3  -  ^31^1 


(41) 


We  let  V^  *  1  and  A|^  ia  the  magnitude  and  phaae  of  the  current  on  element  K. 
Therefore, 

^2  ^21 


Am  “ 


(42) 


Thua,  the  coefficient  A  ^  la  the  aelf  admittance  of  the  active  element  no.  1  and  the 
coefficient  A^.  where  K  la  unequal  to  one,  ia  the  mutual  admittance  between  the 
first  and  K  element. 


2.2  Numerical  Reaulta 

The  simple  case  of  a  =  2X.  and  =  a  +  *./4  (see  Figure  7)  was  investigated 
and  the  mutual  and  self  admittances  were  determined  through  use  of  computer 
programs  to  numerically  solve  the  m  nonhomogeneous  simultaneous  equations 
(Eq.  (37)).  Equation  (37)  becomes 


ni 

gc 

E 

E 

m  =  1 

n=0 

in(dk+b/r^-d^) 


BESS(n)  = 


~V^ 

2PoC 


(43) 
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where 


k  « 

BESSM-J  daQ(o)  Rp¥^(n.a)  +  y  do  Q(o)  ReFB(n, ») 
0  ke« 


+ 


do  Q(o)  lmPS(n,  o)  + 


•e 

/ 


do  Q(o>  lmPB(n,  o) 


} 


(44) 


Q(o) 


(45) 


ReFS(n.o)  •  Jn(x)J„(x) 


J^(x)j;(y)  -  Y;(x)j;(y)  2J„(x)Y^(x)J^(y)Y^(y) 


jS(y>  +  vj(y) 


(46) 


RePB(n.o)  »  yx)Kj^(x) 


(47) 


ImFS(n.a)  «  J„(x)Y„(x) - s - -5 

"  "  Jn<y>-^vy> 


(48) 


2  j  ' 

ImPB(n..)  .  I 


(49) 


For  Eq.  (43)  to  be  solvable,  the  expression 


^  ‘n(4k+*>/r^-4^) 

2^  e  °  BESS(n) 

n=0 


(50) 


must  converge  to  aero  for  large  n,  otherwise,  the  appropriate  coefficients  for  each 
(m  =  1,  2,  . . .  m)  cannot  be  determined.  Let 


SR(n)  =  /Qfe)ReFS(n,o)  do 
SI(n)  *  /Q(o)ImFS(n,o)  4» 
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BR<n)  *  /Q6>)ReFB(n.o)  ob 
Bl(n)  =  /Q(ar)ImFB(n.cr)  (far 

The  col  trlbutions  to  BESS(n)^^  of  SR  and  BI  rapidly  dlmlnlah  for  large  order  n. 

BR  and  SI  decrease  very  slowly  In  absolute  Vt^lue  for  Irrge  order  n.  Figures  8,  9, 
10,  and  11  show  the  values  of  SR,  SI,  BR  and  BI  respectively  for  order  n  =  0  to 
700.  BR  and  SI  approach  sero  very  slowly  for  large  order  n,  while  SR  and  SI  con¬ 
verge  rapidly  to  sero  for  large  n.  Figure  12  shows  the  value  BSlJMln)  where 

BSUM(n)  =>  '^{SR(n)  +  BR(n))*  +  (SJ(n>  +  Bl(n))*  *  iBESSInlj 


and  (worst  case).  BSUM  Is  shown  to  slowly  converge  to  sero  for 

large  order  n.  Due  to  time  limitation  In  computer  usage,  n  ==  0  to  700  was  chosen; 
however,  for  more  reliable  results,  n  >  0  to  2000  may  be  required,  especially  for 
the  worst  case  =  0°. 

Once  l.\c.'  BCSS(n)  were  determined  for  n  =  0  to  700,  these  values  were  used  In 

Eq.  (50)  to  determine  the  coefficients  for  the  A_ 's  in  the  set  of  simultaneous  non- 

m 

homogeneous  equations  (Eq.  (37)).  Each  resultant  matrix  was  found  to  be  of  the 
form 


"^l  '^'m  •  •  •  *^2 

!'*i' 

^1 

Tg  Tj  ...  Tq 

^2 

•  f  *  « 

^3 

•  «  •  « 

^4 

^4 

•  •  «  « 

• 

• 

T  T  ,  V  ,  •  •  T, 

A 

V 

m  ni  - 1  1 

m 

m 

(51) 


Thus,  due  to  symmetry,  Eq.  (37)  need  only  be  evaluated  once  for  any  m  =  l,  2,  ... 

m  to  determine  all  coefficients  for  the  A  's  in  the  equation. 

m 

Results  of  the  self  and  mutual  admittances  for  a  circular  array  of  cylindrical 
dipoles  a  quarter  of  a  wavelength  above  a  perfectly  conducting  cylinder  of  radius 
r  -  2X  were  calculated.  See  Figure  12  for  array  sizes  of  m  =  2  to  22.  Figure  13 
shows  the  self  admittance  for  the  active  element,  Figure  14  shows  the  mutual  ad¬ 
mittance  between  dipole  no.  I  and  dipole  no.  m.  Due  to  symmetry,  these  two 
should  be  the  same.  The  inconsistency  may  be  the  result  of  not  taking  large  enough 


ll.  Amos,  V.  E, ,  and  Daniel,  S,  L,  (1977)  AM03LIB,  A  Special  Library  Version. 
Sandla  Laboratories  (SAND  77-1390). 
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order  n  in  Eq.  (50)  to  have  BESS(n)  approximately  equal  to  zero.  Figure  15  shows 
the  mutual  admittance  between  dipole  no.  1  and  the  adjacent  dipole  no.  m.  Fig¬ 
ure  16  shows  the  relative  currents  for  the  dipoles  for  arrays  of  size  m  =  2,  4,  7, 
10,  13,  16.  Once  again,  for  m  s  12,  asymmetries  occur  which  may  be  due  to  not 
taking  high  enough  order  n  for  the  evaluation  of  BESS(n)  in  Eq.  (50)  for  the  deter¬ 
mination  of  the  coefficients  for  the  A^^'s. 


Figure  7.  Array  of  Cylindrical  Dipoles  About  a  Perfectly  Conducting 
Cylinder,  Simple  Case,  a  =  2X,  p  =  =  a  +  X/4 
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Figure  10.  BR(n)  =  /  Q(ar)  ReFB(n,  o)  for  n  =  0  to  700 


Figure  11.  BI(n)  -  /Q(o)  ImFB(n,  o)  for  n  =  0  to  700 


ELEMENT  NUMBER 


ELEMENT  NUMBER 


Figure  16.  Relative  Currents  For  the  Circular  Array  of 
Dipoles  Above  a  Perfectly  Conducting  Cylinder  For  Active 
Element  Dipole  No.  1,  m  =  2.  4,  7,  10,  13,  16 
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Figure  16.  Relative  Currents  For  the  Circular  Array  of 
Dipoles  Above  a  Perfectly  Conducting  Cylinder  For  Active 
Element  Dipole  No.  1,  m  =  2,  4,  7,  10,  13,  16  (Cont. ) 
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S.  CONCLUSION 


It  haa  been  ahown  that  numerical  techniquea  for  approximation  of  the  aolutlon 
of  the  m  nonhomogeneoua  aimultaneoua  equatlona  (Eq.  (37))  provide  uaefUl  and  ana¬ 
lytical  data  In  the  determination  of  characteristics  of  a  circular  array  above  a 
perfectly  conducting  cylinder.  From  Eqs.  (37)  and  (50),  it  is  seen  that  the  value 
BESS(n)  must  converge  to  xero  for  large  order  n  to  determine  the  coefficients  of 
the  A^'s  in  Eq.  (39).  Due  to  the  symmetry  of  the  coefficients,  as  seen  in  Eq.  (51), 
the  values  T|2<  •  •  •  need  only  be  calculated  once  for  a  given  array  of  size 
m.  Also,  from  the  expression  for  BESS(n),  Eq.  (50),  BESS(n)  is  Independent  of  m, 
the  number  of  dipole  elements  in  the  circular  array  being  considered.  BESS(n)  is, 
however,  dependent  on  the  values  of  the  radius  of  the  cylinder  and  the  frequency 
(Wavelength).  Consequently,  tables  can  be  made  of  the  BESS(n)  values  for  n  =  0 
to  N,  where  N  is  of  order  n  large  enough  to  force  BESS(n)  to  zero.  These  tables 
would  enable  calculation  of  the  mutual  and  self  admittances  for  an  array  of  arbitrary 
size  m  about  a  perfectly  conducting  cylinder  of  radius  r'.  An  understanding  of  the 
effects  of  mutual  coupling,  taking  into  consideration  both  the  space  wave  and  the 
creeping  wave,  is  necessary  in  the  synthesizing  of  an  aperture  distribution  to  give 
the  best  fit  to  a  specified  radiation  pattern,  and  for  the  determination  of  optimiza¬ 
tion  techniques  for  circular  array  antennas. 
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